Time-Optimal Transfer of Coherence 



(N 

o 

(N 

> 
O 



43" 
^— > ■ 



> 
in 
o 

p 

OV 

o: 



X 



Alberto Carlini 1 ' 2 ^ and Tatsuhiko Koikc 3 'Q 

1 Dipartimento di Scienze ed Innovazione Tecnologica, 
Universita' del Piemonte Orientale, Alessandria, Italy 
2 Istituto Nazionale di Fisica Nucleare, Sezione di Torino, Gruppo Collegato di Alessandria, Italy 
'Department of Physics, Keio University, Yokohama, Japan 
(Dated: November 10, 2012) 

We provide exact analytical solutions for the problem of time-optimal transfer of coherence from 
one spin polarization to a three-fold coherence in a trilinear Ising chain with a fixed energy available 
and subject to local controls with a non negligible time cost. The time of transfer is optimal and 
consistent with a previous numerical result obtained assuming instantaneous local controls. 
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I. INTRODUCTION 

The techniques of quantum optimal control are ubiq- 
uitous and have multiple applications in, e.g., molecular 
processes nuclear magnetic resonance (NMR) spec- 
troscopy [J quantum information theory I n par- 
ticular, time-optimal control, where the goal is to min- 
imize the physical time to reach a quantum target (i.e. 
a quantum state or a unitary operation) is important 
for the construction of efficient gates in quantum com- 
puting architectures and it offers a more physical frame- 
work in defining the complexity of quantum algorithms 
[jj. On the other hand, controlling the spin dynamics in 
quantum chains is of relevance in the efficient exchange 
of quantum information [f|, and, e.g., in NMR multi- 
dimensional spectroscopy experiments. The problem of 
the efficient transfer of coherence in a trilinear spin chain 
where neighbor spins are subject to equal Ising couplings 
and where the single qubits are separately addressed via 
instantaneous, local controls was discussed within a ge- 
ometrical quantum control ansatz in Q. The case of 
unequal Ising couplings and to the case of longer chains 
was treated in The time duration of the trans- 
fer given numerically in Q , though shorter than the one 
obtained with conventional approaches, was not guaran- 
teed to be a (global) time-optimal control. One of the 
postulates in the above works is that the time cost of 
one-qubit operations is zero, which is reasonable (e.g., 
in heteronuclear NMR) when the other time-scales (e.g., 
the inverse of the maximal Ising couplings) involved in 
the quantum control process are much longer. In a se- 
ries of papers [IH-[l5|, we discussed a theoretical scheme 
for time-optimal quantum control, which we called the 
quantum brachistochrone (QB [12]). where the postu- 
late of zero-time local operations is not required. The 
scheme was The QB is derived from an action principle 
which enforces the dynamical laws of quantum evolution 
(i.e., the Schrddingcr equation or a master equation) and 
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the constraints which the Hamiltonian of the physical 
system has to satisfy (e.g., a fixed total energy, the ab- 
sence of certain qubit interactions etc.). The framework 
was designed for the time-optimal evolution of quantum 
states between fixed initial and final states [l3j]-[lH, for 
the time-optimal generation of a certain unitary quantum 
gate [l5|, and for the situation (typical in experiments) 
where the target is reachable only in an approximate way 
(i.e., with a fidelity smaller than one (l6l|). 

In this brief note we reconsider the problem of the time- 
optimal transfer of one qubit polarization to three qubit 
quantum coherence studied in Q. We use a trilinear 
qubit system with the same Ising interaction Hamiltonian 
as in [7[, but we consider a generic and time consuming 
local control on the intermediate qubit of the chain and 
we assume that a finite energy is available. By using 
the QB methods, we are able to analytically quantify the 
shape of the local control and the duration of the optimal 
transfer, which essentially coincides with the numerical 
value found in Q. 



II. QUANTUM BRACHISTOCHRONE 

Let us first review the main formalism of the QB as dis- 
cussed, e.g., in [l5|. The problem of the QB for unitary 
operations is to find the time-optimal way to generate a 
target quantum gate Uf via the control of a Hamiltonian 
H(t). We impose that the quantum evolution is driven 
by the Schrddingcr equation, and that the Hamiltonian 
has to satisfy certain constraints. For example, the en- 
ergy available in the experiments may be limited, or some 
interactions between the qubits may be forbidden. The 
QB problem can be concisely formulated in terms of the 
action [H,[l5[: 

S(I7 J fr;a ) A ) A i ) := / ds [Na + L s + L c ] , (1) 
Jo 

L s :=(A,i%tf -aH), (2) 

L c :=aJ2^f j (H), (3) 



koikc@phys.kcio.ac.jp 



2 



where (A,B) := ^(A^B) and logA^ is the number of 
qubits. The Hermitian operator A(s) and the real func- 
tions Xj(s) are Lagrange multipliers. The function a(s) 
gives the time cost and it relates the physical t and the 
parameter s times via t := J a(s)ds [14| . 

The multiplier A enforces the Schrodinger equation: 



dU 
l ~dt 



HU. 



(4) 



The constraints on the Hamiltonian are obtained from 
the variation of the action (TTJ) with respect to Xj : 

fj(H)=0. (5) 

For instance, the finite energy condition reads: 



f (E) :-i[Tr(iJ 2 ) 



Nu 2 



0, 



with uj := const. 

The variation of S with respect to E gives: 



A = X H 



^ x d fj (E) 



^ 1 BE 



(6) 



(7) 



while the variation with respect to a gives a normaliza- 
tion condition for A: 



Tr(ifA) = N. 



(8) 



Finally, variation of the action with respect to U gives, 
after some elementary steps: 



dA 
dt 



[E,A\. 



(9) 



We name eq. (JSJ) the quantum brachistochrone equa- 
tion. Together with the constraints <j5j) it defines a 
boundary-value problem for the evolution operator U(t) 
with fixed initial (U(t = 0) = 1) and final conditions 
(U(t = T) = Uf), where T is the optimal time duration 
of the quantum evolution. Once the target quantum gate 
Uf is given, one can solve the quantum brachistochrone 
© together with the constraints and find the time- 
optimal Hamiltonian H opt (t), which will depend upon 
a set of integration constants. Then, integration of the 
Schrodinger equation d3| from U(0) = 1 produces the 
time-optimal evolution U op t(t). Finally, the integration 
constants (including the time duration T ) in E opt {t) are 
determined imposing that /7 op t (T) — Uf [ 1 51 ] . 



III. ISING HAMILTONIAN 

Let us now explicitly use the QB formalism for the 
model of a linear Ising chain of three qubits (identified 
by indices a G {1,2,3}) with Ising couplings J 12 and 
J23. We further assume that the intermediate qubit of 
the chain is controlled via a local magnetic field Bi(t) 
(i = x,y,z), i.e., we consider the Hamiltonian: 



E(t) := Ji 2 ala 2 z + J^a^o* + B(t) ■ a 2 , 



(10) 



where, e.g., o\a 2 := cr, (8) Oj ® 1, and CTj are the Pauli 



operators [3JJ. The interaction part in eq. (fTUj) is the 
same as E c in eq. (2) of Q. However, the local and 
time consuming term B(t) ■ a 2 is used instead of the zero 
cost terms Ha and Hb of 0- We also introduce the 
couplings' ratio K := J23/J12, rescale time as r := J\$b, 
the magnetic field as B{t) := B(t)/J 12 and the energy 
as u) := u)j J\2- Then, e.g., the finite energy condition © 
explicitly reads: 



B = lu 2 - (1 



K 2 ) = const. 



(11) 



The analytical solution of the quantum brachis- 
tochrone problem for the Hamiltonian (fTU)) was found 
in Ref. [16| (for details, se Sections III and VIIIA there). 
Here we are interested in one of the main results, i.e. the 
time-optimal magnetic field which reads: 



Sopt(r) 



Bo cos 9(t) 
B sin 6(t) 
B z 



(12) 



where 8(t) := (It + 9q and fl, B z , Bq and are integra- 
tion constants. 

The explicit form of the time-optimal evolution opera- 
tor U op t(t), which in principle can be obtained integrating 
the Schrodinger eq. (|4]) and using eqs. (fTU)) and ([12]) . is 
not necessary for the purposes of this paper (for details 
see [13). 



IV. TIME-OPTIMAL TRANSFER OF 
COHERENCE 

Our goal now is to time-optimally create a 3-qubit 
quantum coherence starting from an initial polarization 
on one of the qubits at the ends of the Ising chain. For 
instance, we would like to realize the following transfor- 
mation: 



123 



(13) 



and find the time-optimal control of the Hamiltonian ([10]) 
which performs this transformation. We call this as the 
time-optimal coherence transfer (TOCT) problem. 

In the previous literature 0, the TOCT problem 
has been addressed in the following way. First, one 
introduces expectation values of operators (0)(t) := 
Tr[Op(r)], where the density operator p(r) = 
U(t)p(0)U'(t) and U(t) satisfies the Schrodinger equa- 
tion ([3]) with the Hamiltonian (fT0"|) . Then, one notes that 
the algebra generated by the Hamiltonian (fTU|) and the 
operator a\ is just an 8-d subspace of the full 63-d algebra 
for the three qubits, and is spanned by x\ := (cr^), X2 



x 6 := (a},ala 3 



. 'XI 5 

{(Tltf), x 5 := (crla 3 



zn 



lT V <T z (J z) ano - ^8 

i([H,0]) one gets a system of 



Using (0)- = Tr[Op] 
equations for the XiS, which depend upon the unknown 
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control B(t) in (|10p . Then one should solve such a sys- 
tem and determine the optimal control B op t{r) ensuring 
that the transfer from the initial condition x\{Q) = 1 to 
the final condition x%(T t ) = 1 happens in the minimal 
time r*. This method has been applied in as- 
suming that one can apply a sequence of selective hard 
pulses on individual qubits in the chain, and numerical 
solutions have been found for the time duration t* for an 
arbitrary couplings' ratio K. 

Here, however, we proceed in a different way. In fact, 
from the variational principle and the action ([!}, we al- 
ready have proof that the general time-optimal unitary 
evolution (up to an arbitrary target) must satisfy the QB 
equation §§§ , and that such evolution is generated by the 
Hamiltonian (flU)) with optimal magnetic field (TT^j) . The 
global time-optimal solution depends upon the specific 
boundary conditions and can be found by consequently 
determining the integration constants in (|12[) . Therefore, 
we directly insert KTUJ and (H2J) into (O)' =i([H,0}) and, 
defining the vector x := (x\, xi, 23, 24, X5, x§, xj, xg)*, we 
obtain the system of equations: 



dx{r) 
dr 



where the 8x8 matrix: 



M(r) 



M(r)x(r), 



P(r) Q 
Q P(r) 



(14) 



(15) 



depends on the skew symmetric, time-dependent 4x4 
matrix: 



P(t) :- 



- 





-1 











Bq sin 9 


-B z 


1 


—Bq sin 9 





Bq COS 


. 


B z 


— Bq cos 9 






(16) 



(with 9(t) given below eq. (fT^j) ) and on the skew sym- 
metric, constant 4x4 matrix: 



Q := K 




-1 

10 



(17) 



Our task is to solve the system (jT4")) subject to the bound- 
ary conditions xi(0) = 1,Xs(t*) = 1, and to deter- 
mine the values of the unknown integrals of the mo- 
tion Bo , B z and Cl in such that the time duration 
of the transfer t* is minimal. We simplify the prob- 
lem by introducing the 4-dimensional vectors x+(t) := 
{xi,X2,x%,X4) t and x_(t) := (x 5 , x§, 27, Xg)*, their com- 
binations: y±(r) :~ x + (t)±X-(t) and the 4 x 4 matrices 
M ± (t) := 2[P(t)±Q}. Then, cq. (TJH) decouples into the 
two sets of first order, linear differential equations: 



dy± (t) 
dr 



M±(r)y ± (r), 



(18) 



whose solutions are: 

V±(r) = cxp[A ± (r)]y±(0), 
where the 4x4 matrix: 



(19) 



A ± (t) 
and the 2x2 matrices 



R±(t) 



^ [cos 6»(t) -cos6»o] (B z ±K)t 



Bn 



Rq(t) := i-^-[sin0(r) — sin9o]a y . 



, (21) 



(22) 



For the transfer of coherence (|13p . the initial and fi- 
nal boundary conditions can be rewritten as y+{0) = 
£1(0) := '(1,0,0,0)* and y+{n) = -yl(r*) := 
(0, 0, 0, 1)'. After a long but straightforward analysis (for 
a sketch of the derivation see the Appendix), we finally 
find that the time-optimal transfer of coherence is char- 
acterized by the following integrals of the motion. The 
minimal time duration of the transfer is II 811 : 



T*opt 



V3 



-7T, 



(23) 



while the optimal magnetic field (also using the energy 
constraint (j5])) reads: 



B z 
Bo 
O 







±-\[gF- 

TT 



9o = ^ [(2r + 1)tt ± \/3 V^ 2 - 2] 



(24) 
(25) 
(26) 

(27) 



where r is an arbitrary integer and Cj > 2. It is worth 
noticing that the same time-optimal solution (with time 
duration and control as in (|23[) -([27 |) ) is also valid for for 
the transfer of coherence from the initial polarization a\ 
to the final operator (jyO^o^. [lj|. No transfer is possible, 
instead, from <j\ to (jyO^a?,. 



V. DISCUSSION 

In this paper we described the most recent develop- 
ments of the QB formalism. We investigated the model 
of a trilincar Ising Hamiltonian with equal (or opposite) 
Ising couplings, subject to a local, time consuming con- 
trol on the intermediate qubit and with a finite energy 
available. The formalism of the QB [l3| enables to find 
the exact, analytical form of the laws of the time-optimal 
control for the transfer from a single spin polarization a\ 
to the final three-qubit coherence a) g a^ j a\. The analytical 
expression of the time duration of the transfer, equation 
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(l23|) . matches with a numerical result which appeared in 
the previous literature 0] for the same Ising interactions 
in the trilinear chain, but with only instantaneous local 
controls available (20[. It was proven here that this time 
of transfer is optimal. An extension to the case of differ- 
ent couplings in the trilinear chain, to the situaiton where 
local controls are allowed on all the qubits in the chains, 
and to chains with more qubits and nonlinear topologies 
is in progress. It is also important to consider the more 
general ansatz where the total energy available is not a 
constant, but only bounded from above. 

In the standard paradigm of time-optimal quantum 
computing (see, e.g., fTT| ) one-qubit unitary operations 
are assumed to have zero time cost. In more physical 
situations where this time cost cannot be neglected with 
respect to the other time-scales appearing in the quan- 
tum control problem, the neat geometrical methods of 
the standard paradigm cannot be immediately applied. 
The QB provides a natural formalism for this scenario. 



and introducing the quantities 
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APPENDIX 

The initial and final boundary conditions y+(0) = 
y_(0) := (1,0,0,0)* andy+(n) = -y_(Y*) := (0,0,0,1)* 
can be rewritten in terms of the matrix A± as: 



,A±(r„) 



31 



±1. 



41 



(28) 



Using the notation: 



a := 2t„, 
c± := (B z ± K)a, 

b := 2^^[cos0(r*) — cos(9 ], 
d:=2^[sin0(r*) -sin0 o ], 



(29) 
(30) 

(31) 
(32) 



(33) 
(34) 

(35) 

(36) 
(37) 

(38) 



we can explicitly write the boundary conditions (|28[) (in 
the order of components (11)±; (31)±; (21)±; (41)±) as: 





A+]C| = [X+ - 2a 2 
A~]CY = [X- - 2a 2 



<A + )Si = [X+ - 2< - 
= [X- - 2c 2 . - 
b(Cf - C z ) = c + d(Sf - S z ), 




b(C™ - C™) = c-d(S™ - 55), 



(40) 

'A+]Sf (41) 
^AT]^ (42) 
(43) 
(44) 



d{C + 



' C z ) 



-c+b(Sf - S z ) + 



-c-b(Sl 



S 



(45) 
(46) 



An elementary analysis of the system of equations (f3"9")l - 
(|4^]) leads to find out that its only solution is possible for 
the values of the parameters d — and b ^ 0. We now 
sketch how we get to the complete, analytical form of 
the solutions for the time-optimal transfer of coherence. 
First, for d = but 6^0, and from eqs. one 
finds that Z_ = — Z + +2np and W + = —W-+2itq, where 
p and q arc arbitrary and (du e to the definitions (|35[1 - 



nonncgativc integers (2 lj . Then, substituting the 



latter relations into (|39p - P0|) . we find that the only pos- 
sible choice H| is Z+ = f (2n+l) and W- = §(2m+l), 
where again m ad n are arbitrary, nonnegative integers. 

We can now eliminate Z± and W± from the square of 
eqs. ([33 ]) -([55 1) to get: 




2ir 2 p[(2n + 1) - 2p], 
-2n 2 q[(2m + l) - 2q], 



(47) 
(48) 



-[(2n + l) 2 -4p(2n + l) + 8p 2 ], (49) 



a 2 + b 2 + c 2 _ = —[(2m + l) 2 - 4g(2m + 1) + 8q 2 ]. 



(50) 



In other words, we have replaced the original system 
of eqs. ([3"9" l) -([4l) |) with eqs. (|^7 ]) -([5D ]) (instead of eqs. 
-([iD]) and eqs. (|S5 ]> - (flS jl ) and (after substitution of 
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Z_,W_ into cqs. (|3I ]) -([i2 ]l and (|4*S" ]) -([4"6 ]) ) with equa- After a short algebra, the solution to eqs. d47j) - d54j) is 
tions: given by 2p = 2q = ni + n+ 1 and the following constants: 

7A7 = (a 2 + 6 2 - 4) „„ -. , (51) 



[(2n + 1) - 2jj] ' a = ±c± = - yf(2m + l)(2n + 1), (55) 



y/Al=-(a 2 + b 2 -c 2 _) J — -^1— ; -, (52) 6=-7rA-(n-m), (56) 

[(2m + 1) — 2q\ 

abc + = (— l) n+1 (a 2 + b 2 - g 2 ^ ^A^+jOK^+jO 4p] ^ |^ _ ^ where the integers to and n have op- 



posite parity (m = 2too and n = 2uq + 1) and n > n 
(i.e. no > toq > 0). Finally, we can use the results (|55|) 



4[(2n + 1) - 2p] 

(53) 

^ c _ / ^m/ a 2 _|_ ^2 c 2\ tt(2to + 1)[(2to + 1) — 4q] (|5rJ|) to invert (j2"5j) - (j3"2")) in favor of the physical quantities 

4[(2to + 1) — 2q] r*, Bo, B z , 0, and 9q. Minimization of t* with respect to 

(54) to and n gives eqs. (f!?3"|) - (|27[) of the main text. 
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